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Abstract. The paper presents an integrated view of the population structure and its role in 
establishing the ionisation state of light eleinents in dynamic, finite density, laboratory and 
astrophysical plasmas. There are four main issues, the generalised coUisional-radiative picture 
for metastables in dynamic plasmas with Maxwellian free electrons and its particularising 
to light elements, the methods of bundling and projection for manipulating the population 
equations, the systematic production/use of state selective fundamental collision data in the 
metastable resolved picture to all levels for collisonal-radiative modelling and the delivery of 
appropriate derived coefficients for experiment analysis. The ions of carbon, oxygen and neon 
are used in illustration. The practical implementation of the methods described here is part of 
the ADAS Project. 
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1. Introduction 

The broad mechanism for radiation emission from a hot tenuous plasma is simple. Thermal 
kinetic energy of free electrons in the plasma is transferred by collisions to the internal energy 
of impurity ions, 

A + e-^A*+e (1) 

where A* denotes an excited state and A the ground state of the impurity ion. This energy is 
then radiated as spectrum line photons which escape from the plasma volume 

A* ^A + hv. (2) 

where hv is the emitted photon energy and v its frequency. Similarly ions in general increase 
or decrease their charge state by colhsions with electrons 

A + e A^ + e + e 

A" +e^A + hv (3) 

where A^ denotes the next ionisation stage of impurity ion A. The situation is often referred 
to as the coronal picture. The coronal picture has been the basis for the description of 
impurities in fusion plasmas for many years. However, the progress towards ignition of 
fusion plasmas and to higher density plasmas requires a description beyond the coronal 
approximation. Models of finite density plasmas which include some parts of the competition 
between radiative and collisional processes are loosely called colUsional-radiative. However, 
coUisional-radiative theory in its origins (Bates et al , 1962) was designed for the description 
of dynamic plasmas and this aspect is essential for the present situations of divertors, heavy 
species, transport barriers and transient events. The present work is centred on generalised 
collisional radiative (SC3^) theory (McWhirter and Summers, 1984) which is developed in the 
following sections. It is shown that consideration of relaxation time-scales, metastable states, 
secondary collisions etc. - aspects rigorously specified in collisional-radiative theory - allow 
an atomic description suitable for modelhng the newer areas above. The detailed quantitative 
description is complicated because of the need to evaluate individually the many controlUng 
collisional and radiative processes, a task which is compounded by the variety of atoms and 
ions which participate. The focus is restricted to plasmas which are optically thin and not 
influenced by external radiation fields, and for which ground and metastable populations of 
ions dominate other excited ion populations. The paper provides an overview of key methods 
used to expedite this for light elements and draws illustrative results from the ions of carbon, 
oxygen and neon. The paper is intended as the first of a series of papers on the application of 
collisional-radiative modeUing in more advanced plasma scenarios and to specific important 
species. 

The practical implementation of the methods described here is part of the ADAS (Atomic 
Data and Analysis Structure) Project (Summers, 1993, 2004). Illustrations are drawn from 
ADAS codes and the ADAS fundamental and derived databases. 

1.1. Time constants 

The lifetimes of the various states of atoms, ions and electrons in a plasma to radiative 
or collisional processes vary enormously. Of particular concern for spectroscopic studies 
of dynamic finite density plasmas are those of translational states of free electrons, atoms 
and ions and internal excited states (including states of ionisation) of atoms and ions. 
These hfetimes determine the relaxation times of the various populations, the rank order of 
which, together with their values relative to observation times and plasma development times 
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determines the modelling approach. The key lifetimes divide into two groups. The first is 
the intrinsic group, comprising purely atomic parameters, and includes metastable radiative 
decay, t^, ordinary excited state radiative decay t„ and auto-ionising state decay (radiative 
and Auger), Ta- The intrinsic group for a particular ion is generally ordered as 

Ta « To « T,„ (4) 

with typical values 

T,„ ~ 10 S,To ~ lO^Vz"* S,Ta ~ 10"'^ S. (5) 

where z is the ion charge. The second is the extrinsic group, which depends on plasma 
conditions - especially particle density. It includes free particle thermalisation (including 
electron-electron r^-e, ion-ion t,_,- and ion-electron t,_(.), charge-state change (ionisation t,o„ 
and recombination t,.j,c) and redistribution of population amongst excited ion states (Tred)- The 
extrinsic group is ordered in general as 

Tion.rec » » T,_; >> T^-e (6) 

with approximate expressions for the time constants given by 

Tree ~ [10" - 10''](1/(Z+ l)2)(^r,/fe)'^'(cm-VNe) S 

Tio„ - [10' - IQ'Kz + l)4(/HAr,)i/V^*^'(cm-VNe) s 
T,_,- ~[7.0xl0'](m,/m^)'^2 

(kTe/lH)"\l/z'')(cm-^m s (7) 
r,_, ~[1.4xl0^](m,/m^)i/2 

((kTe/In) + 5.4 X lO-\kTi/lH)(mp/mdf'^ 
(l/z2)(cm-VNi) s 

Te-e ~ [1.6xlO''](kTe/lH?'\cm-'m S. 

The ion mass is m,, the proton mass nip, the ionisation potential the ion density A^,, the 
electron density A^^., the ion temperature J,, the electron temperature and the ionisation 
energy of hydrogen is In- Tred may span across the inequalities of equation|6]and is discussed 
in a later paragraph. 

From a dynamic point of view, the intrinsic and extrinsic groups are to be compared with 
each other and with timescales, Tpiasma, representing plasma ion diffusion across temperature 
or density scale lengths, relaxation times of transient phenomena and observation times. For 
most plasmas in magnetic confinement fusion and astrophysics 

where Tg represents the relaxation time of ground state populations of ions (a composite of 
Tree and T,o„) and it is such plasmas which are addressed in this paper. These time-scales 
imply that the dominant populations of impurities in the plasma are those of the ground and 
metastable states of the various ions. The dominant populations evolve on time-scales of 
the order of plasma diffusion time-scales and so should be modelled dynamically, that is in 
the time-dependent, spatially varying, particle number continuity equations, along with the 
momentum and energy equations of plasma transport theory. Illustrative results are shown in 
figure [^. 

The excited populations of impurities and the free electrons on the other hand may 
be assumed relaxed with respect to the instantaneous dominant populations, that is they 
are in a quasi-equilibrium. The quasi-equilibrium is determined by local conditions of 
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Figure 1. (a) Ratios of ground and metastable lifetimes to plasma timescales arising from 
transport of some ions of nitrogen and neon at the edge of fusion plasmas. For N" the ground 
state to metastable state transfer time constant, t„ is contrasted with t,, = A,„i/v with typical 
scrape-off-layer thickness /i,,,/ ~ 2 cm and mv^ /2 < 1 eV from chemical or physical sputtering 
and with t^,/ = /i^^j/f where a typical diffusion coefficient D = lO** cm^s"'. For selected 
neon ions, tj = ir^.w, /I'rfi/ is contrasted with the reciprocal of the sum of the inverses of 
the ionisation and recombination time constants for the ground state. The latter illustrative 
results are for a JET-like tokamak H-mode radial plasma model with shaped diffusion and 
pinch terms. T^ir = 0)=1 keV with a pedestal of 30eV at r = a and exponential decay (scale 
length = O.Ola) in the scrape of layer at plasma minor radius a = 100 cm. iv,/ combines 
pinch and concentration diffusion parts, (b) Critical densities, defined as NgiTred = tq.o), for 
categories of excited populations of some oxygen ions. The x-axis scale for curve [1] is the 
orbital angular momentum / of the outer electron of doubly excited states; for curve [2] it is the 
principal quantum number n of the singly excited electron; for curve [3] it is a simple index 
to the three low-lying states illustrated. For low levels of ions, T,-fd is markedly sensitive to 
the detailed atomic structure, n = 3 valence shells populations are of special relevance to light 
element spectroscopy in the visible. 
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electron temperature and electron density. So, the atomic modelling may be partially de- 
coupled from the impurity transport problem into local calculations which provide quasi- 
equilibrium excited ion populations and emissivities and then effective source coefficients 
(collisional-radiative coefficients) for dominant populations which must be entered into the 
plasma transport equations. The solution of the transport equations establishes the spatial and 
temporal behaviour of the dominant populations which may then be re-associated with the 
local emissivity calculations for matching to and analysis of observations. 

For excited populations, r^j plays a special and complicated role due to the very large 
variation in collisional excitation/de-excitation reaction rates with the quantum numbers of 
the participating states. In the low density coronal picture Tred » Tq and redistribution plays 
no part. Critical densities occur for Tred ~ and for Tred ~ and allow division of the 
(in principle) infinite number of excited populations into categories including low levels, high 
singly excited levels and doubly excited levels for which important simplifications are possible. 
These are examined in section |2 Light element ions in fusion plasmas are generally in the 
singly excited state redistibutive case, approaching the doubly excited redistribution case at 
the higher densities. Highly ionised ions of heavy species in fusion plasmas approach the 
coronal picture. Illustrative results on critical densities are shown in figuren^- 

Finally, because of the generally short Te-e compared with other timescales (including 
those of free-free and free-bound emission), it is usually the case that the free electrons have 
close to a Maxwellian distribution. This assumption is made throughout the present paper, but 
is relaxed in the next paper of the series (Bryans et al , 2005). 

1.2. Generalised collisional-radiative theory 

The basic model was established by Bates et al (1962). The ion in a plasma is viewed 
as composed of a complete set of levels indexed by / and j and a set of radiative and 
collisional couplings between them denoted by C,j (an element of the collisional-radiative 
matrix representing transition from j to /) to which are added direct ionisations from each 
level of the ion to the next ionisation stage (coefficient S,) and direct recombinations to each 
level of the ion from the next ionisation stage (coefficient r,). Thus, for each level, there is a 
total loss rate coefficient for its population number density, A',-, given by 



Following the discussion in the introduction, it is noted that populated metastable states can 
exist and there is no real distinction between them and ground states. We use the term 
metastables to denote both ground and metastables states. Metastables are the dominant 
populations and so only recombination events which start with a metastable as a collision 
partner matter. We condider the population structure of the z-times ionised ion, called the 
recombined or child ion. The (z -H l)-times ionised ion is called the recombining or parent ion 
and the (z - l)-times ionised ion is called the grandchild. The metastables of the recombined 
ion are indexed by p and cr, those of the recombining ion by v and v' and those of the 
grandchild by yu and jj.' . Therefore the ion of charge state z has metastable populations Np, 
the recombining ion of charge (z -i- 1) has metastable populations N^, and the grandchild ion 
of charge (z - 1) has metastable populations N'^^. We designate the remaining excited states 
of the z-times ionised ion, with the metastables separated, as ordinary levels for which we 
reserve the indices / and j and populations A^, and A^^. There are then, for example, direct 
recombination coefficients r,-,v from each parent metastable into each child ordinary level and 
direct ionisation coefficients from each child ordinary level to each parent metastable Sv,i such 




(9) 
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that Si — YjvSv,i- Also there are direct ionisation coefficients Sp^^'to the metastables of the 
child from the metastables of the grandchild. Then the continuity equations for population 
number densities are 
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where the equations for the (z - l)-times and (z + l)-times ionised ions have been simplified 
by incorporating their ordinary population contributions in their metastable contributions 
(shown as script capital symbols) as the immediate focus is on the z-times ionised ion. This 
incorporation procedure is shown explicitly in the following equations for the z-times ionised 
ion through to equations ^] and may be done for each ionisation stage separately. Note 
additionally the assumption (made by omission of the (3,1) partition element, where 3 denotes 
the row and 1 the column) that state-selective ionisation from the stage (z- 1) takes place only 
into the metastable manifold of the stage z. 



1.2.1. Derived source term coefficients From the quasi-static assumption, we set dNi/dt - 
and then the matrix equation for the ordinary levels of the z-times ionised ion gives 

Nj = - C^' Ci^N^ - NeCji' riy,N^, (11) 
where we have used summation convention on repeated indices. Substitution in equations II 01 
for the metastables of the z-times ionised ion gives 

dNp r T 

-^^n,[Sp,]n;. 

+ Ne[rpy>-CpjC],'riy,]N^,. (12) 

The left-hand- side is interpreted as a total derivative with time-dependent and convective parts 
and the right-hand-side comprises the source terms. The terms in square brackets in equations 
[21 give the effective growth rates of each metastable population of the z-times ionised ion 
driven by excitation (or de-excitation) from other metastables of the z-times ionised ion, 
by ionisation to the (z + l)-times ionised ion and excitation to other metastables of the z- 
times ionised ion (a negatively signed growth) and by recombination from the metastables 
of the (z + l)-times ionised ion. These are called the SP3? coefficients. Following Burgess 
and Summers (1969), who used the name 'collisional-dielectronic' for 'collisional-radiative' 
when dielectronic recombination is active, we use the nomenclature AGTi for the SC3? 
recombination coefficients which become 

J\.Q*J)y^p = 3^py Tpy CpjCji T/y . (13) 

The SC3? metastable cross-coupling coefficients (for p cr) are 

QCB^^p = tp^lNe = [Cp,r - CpjC^fld^] INe. (14) 

Note that the on-diagonal element \Cpcr - CpjC^j^Cia^ /Ne with cr = p is a total loss rate 
coefficient from the metastable p. Substitution from equation ^2 in ths equations QSI for the 
metastables of the z + 1 -times ionised ion gives 

N,[S„r-SyjC]i'Ci^]N^ 



dt 



+ [Cyy, -N^SyjC]iWiy]N:,. 



(15) 
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The SCIK ionisation coefficients resolved by initial and final metastable state are 

SCDo-^v = §vo- = ['Svo- ~ v/Cj,'C,o-j (16) 

and note that there is contribution to cross-coupling between parents via recombination to 
excited states of the z-times ionised ion followed by re-ionisation to a different metastable, 

■XQ'D,,^,^-Ne[SyjC-jlri,]. (17) 

Consider the sub-matrix comprising the (2,2), (2,3), (3,2) and (3,3) partitions of 
equations[roi Introduce the inverse of this sub-matrix as 



















c 





and note that the inverse of the (1,1) partition 

This compact representation illustrates, that the imposition of the quasi-static assumption 
leading to elimination of the ordinary level populations in favour of the metastable 
populations, may be viewed as a condensation in which the influence of the ordinary levels is 
projected onto the metastable levels. The metastables can be condensed in a similar manner 
onto the ground restoring the original (ground states only) collisional-radiative picture. The 
additive character of the direct metastable couplings Cpa means that these elements may be 
adjusted retrospectively after the main condensations. 



7.2.2. Derived emission and power coefficients There are two kinds of derived coefficients 
associated with individual spectrum line emission in common use in fusion plasma diagnosis. 
These are photon emissivity coefficients CPSC) and ionisation per photon ratio (SXS). The 
reciprocals of the latter are also known as photon efficiencies. From equations the 
emissivity in the spectrum line j —> k may be written as 

M 

ej^k = = Aj^,N,Nj = Aj^kiJ] r'fJ'-'^NeN^ + J] ''^^^^ )■ (20) 

This allows specification of the excitation photon emissivity coefficient 

?£e'"."' ^A/^^r*""^ (21) 

o;j—>k J^'^-' JO- ^ 

and the recombination photon emissivity coefficient 

a'£et::;^, = A;^,^7>. (22) 

The ionisation per photon ratios are most meaningful for the excitation part of the emissivity 
and are 

M,-l 

sxn^^ = £ seD^^jAj^,rj'j''\ (23) 

v=l 

Each of these coefficients is associated with a particular metastable cr, v' or yu' of the A^', 
^+c+i orA^-"' ions respectively. 

The radiated power in a similar manner separates into parts driven by excitation and by 
recombination as 

yL7^ = ^ AEj^kAj^,rj'^''^ (24) 
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called the low-level line power coefficient and 

T3i3v' - 2 ^Ej^kAj^kTf;'^ (25) 
i.k 

called the recombination-bremsstrahlung power coefficient where it is convenient to include 
bremsstrahlung with yXB. Note that in the generalised picture, additional power for the 
z-times ionised ions occurs in forbidden transitions between metastables as 

AE^^pA^^pN^ (26) 

p,o- 

for the z-times ionised ion. In the fusion context, this is usually small. Radiated power is the 
most relevant quantity for experimental detection. For modelling, it is the electron energy loss 
function which enters the fluid energy equation. The contribution to the total electron energy 
loss rate for the z-times ionised ion associated with ionisation and recombination from the 
(z + l)-times ionised ion is 

2 EpiYj ^PcrNcr + NeY^ V^v) + 

p cr v' 

2 AE^^pA^^pN^ + 2 V^^JcrNeN^ + 

p,cr cr 

Y 'P^'ByNeNy (27) 
v' 

where the Ep and Ey are absolute energies of the metastables p of the z-times ionised ion and 
of the metastables v of the (z + l)-times ionised ion respectively. Thus the electron energy loss 
is a derived quantity from the radiative power coefficients and the other generalised coUisonal- 
radiative coefficients. Note that cancellations in the summations cause the reduction to relative 
energies and that in ionisation equilibrium, the electron energy loss equals the radiative power 
loss. 

In the following sections, it is shown how the quasi-static assumption and population 
categorisation allow us to solve the infinite level population structure of each ion in a 
manageable and efficient way. Such solution is necessary for low and medium density 
astrophysical and magnetic confinement fusion plasmas. This is unlike the simpler situation 
of very dense plasmas where heavy level truncation is used, because of continuum merging. 

2. Excited population structure 

The handling of metastables in a generalised collisonal-radiative framework requires a 
detailed specific classification of level structure compatible with both recombining and 
recombined ions. For light element ions, Russell-Saunders (L-S) coupling is appropriate and 
it is sufficient to consider only terms, since although fine structure energy separations may be 
required for high resolution spectroscopy, relative populations of levels of a term are close 
to statistical. So the parent ion metastables are of the form jyEySy with jy the configuration 
and the recombined ion metastables are of the form YpLpSp with the configuration Tp = 
jy + njlj and the excited (including highly excited) terms are {jyLyS y)niliLiS i. n and / 
denote individual electron principal quantum number and orbital angular momentum, L and 
S denote total orbital and total spin angular momenta of the multi-electron ion respectively 
in the specification of an ion state in RusseU-Saunders coupling. The configuration specifies 
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the orbital occupancies of the ion state. For ions of heavy elements, relative populations 
of fine structure levels can differ markedly from statistical and it is necessary to work in 
intermediate coupling with parent metastables of the form jyjy and recombined metastables 
of the form TpJp with the configuration Yp - jy+niU and the excited (including highly excited) 
levels {yyJy)niliiiJi. There is a problem. To cope with the very many principal quantum 
shells participating in the calculations of coUisional-dielectronic coefficients at finite density 
necessitates a grosser viewpoint (in which populations are bundled), whereas for modelling 
detailed spectral line emission, the finer viewpoint (in which populations are fully resolved) is 
required. In practice, each ion tends to have a limited set of low levels principally responsible 
for the dominant spectrum line power emission for which a bundled approach is too imprecise, 
that is, averaged energies, oscillator strengths and collision strengths do not provide a good 
representation. Note also that key parent transitions for dielectronic recombination span a few 
(generally the same) low levels for which precise atomic data are necessary. In the recombined 
ion, parentage gives approximate quantum numbers, that is, levels of the same n (and /) divide 
into those based on different parents. Lifetimes of levels of the same n but different parents 
can vary strongly (for example through secondary autoionisation). Also the recombination 
population of such levels is generally from the parent with which they are classified. We 
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Figure 2. Schematic of population modelling and condensation procedures. It represents two 
metastable parent states of non-zero spin, so that when an nl electron is added, each yields two 
spin systems for the z-times ionised ion. 



therefore recognise three sets of non-exclusive levels of the recombined ion: 

(i) Metastable levels - indexed by p, cr. 

(ii) Low levels - indexed by /, j in a resolved coupling scheme, being the complete set of 
levels of a principal quantum shell range n : uq < n < n\, including relevant metastables 
and spanning transitions contributing substantially to radiative power or of interest for 
specific observations. 

(iii) Bundled levels - segregated according to the parent metastable upon which they are built 
and possibly also by spin system - which can include bundle-nl and bundle-n. 

Viewed as a recombining ion, the set (i) must include relevant parents and set (ii) must span 
transitions which are dielectronic parent transitions. Time dependence matters only for the 
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populations of (i), high precision matters only for groups (i) and (ii) and special very many 
level handling techniques matter only for group (iii). 

To satisfy the various requirements and to allow linking of population sets at different 
resolutions, a series of manipulations on the collisional-radiative matrices are performed 
(Summers and Hooper, 1983). To illustrate this, suppose there is a single parent metastable 
state. Consider the collisional-radiative matrix for the recombined ion and the right hand side 
in the bundle-n picture, and a partition of the populations as [n, n] with n : no < n < n\ 
and «:«!<«. Elimination of the A^,-, yields a set of equations for the N„. We call this a 
'condensation' of the whole set of populations onto the n populations. The coefficients are the 
effective ionisation coefficients from the n, the effective cross-coupling coefficients between 
the n and the effective recombination coefficients into the n, which now include direct parts 
and indirect parts through the levels h. Exclusion of the direct terms prior to the manipulations 
yields only the the indirect parts. Call these C'"if'' and rj"/" . We make the assumption that 
C™J"' and r'^^"' may be expanded over the resolved low level set (see section to give the 
expanded indirect matrix C!"''"' and r!"''"^ where / and j span the resolved low level set (ij). 
These indirect couplings are then combined with higher precision direct couplings Cj"' and 
r^'^ so that 



Cij = Cf/ + qf"' (28) 

and 

'■/V = >f + '•If"- (29) 

The procedure is shown schematically in figure|2] The process may be continued, condensing 
the low level set onto the metastable set. The generalised coUisional-dielectronic coefficients 
are the result. The time dependent and/or spatial non-equilibrium transport equations which 
describe the evolution of the ground and metastable populations of ions in a plasma use these 
generalised coefficients. Following solution, the condensations can be reversed to recover the 
complete set of excited populations and hence any required spectral emission. The progressive 
condensation described above can be viewed as simply one of a number of possible paths 
which might be preferred because of special physical conditions or observations. 

For light element ions, four types of bundling and condensation are distinguished in this 
work: 

(a) Ground parent, spin summed bundle-n — > lowest n-shell. 

(b) Parent and spin separated bundle-n lowest spin system n-shell (the bundle-nS 
population model). 

(c) Low LS resolved metastable states (the low-level population model). 

(d) Parent and spin separated bundle-n — > low LS resolved — > metastable states. 

Type (a) corresponds to the approach used in Summers (1974). Type (c) corresponds to the 
usual population calculation for low levels in which (consistent) recombination and ionisation 
involving excited states are ignored. It establishes the dependence of each population on 
excitation for the various metastables only. Type (d), effectively the merging of (b) and (c) in 
the manner described earlier, is the principal procedure to be exploited in this work for first 
quaUty studies. Details are in the following sub-sections. 
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2s^2p^ 4^ 


1 


N-like 


1 


(2i22;?3 4^-)3„ 


0.375 


2/2/ V 


1 






(2.922753 5„ 


0.625 


2/2/3*55 


1 




2 


(2s^2p^ ^D) 


0.750 


2/2/ ^P 


1 






(2s^2p^ ^D) 'n 


0.250 


2/2/ 'D 


2 




3 


{2s^2p'^ 2p) 3„ 


0.750 


2/2/ 3p 


1 






(2i22;73 2p^ I„ 


0.250 


2/2p4 


2 


O-Uke 


1 


(2^22/ 2„ 


0.250 


2/2p5 2p 


1 






(2i22p4 3f ) 4„ 


0.250 


2*22/35 4p 


1 




2 


(2^22/ !£)) 2„ 


1.000 


2*22^43^ 2£, 


2 




3 


(2*22/ '5) 2« 


1.000 


2/2/3s2£) 


2 


F-like 


1 


(2s22;?5 2p) 1„ 


0.250 


2/2/ 


1 






(2s22p5 2p) 3„ 


0.750 


2s^2p^3s ^P 


1 



Table 1. Bundle-jiS calculation pathways. The parent/spin system weight factor is defined in 
equation l?^ and l?H N,„f,, indicates the number of metastables of the recombined ion associated 
with the parenl/spin system. The parent index, sequentially numbering the different parents 
shown in brackets, is used as the reference for tabulation of coefficients. 



2.1. The bundle-nS model 

Now let yi^'' denote the recombining ion and yi^""'"' the recombined ion so that z = zi - 1 
is the ion charge of the latter. z\ is the effective ion charge and takes the place of the nuclear 
charge in the reduction of hydrogenic rate coefficients to compact forms in the statistical 
balance equations. Also introduce bundled populations 

Nv,nS = N(y^L,S ,),nS - '^N(y^,L^,s ,.),nlLS (30) 
l,L 
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and the assumption that 

(2L+1) _ (2L+1) 

(21 + l)(2Ly + 1) n^(2Lv + 1) 

The bundling is based on the observation that the largest collision cross-sections are those for 
which n - n' and I - I' + \. For these cases the transition energy is small (effectively zero 
for hydrogen or hydrogenic ions) and the cross-sections are so large for electron densities of 
relevance for fusion that it is very good approximation to assume relative statistical population 
for the IL sub-states. The assumption is weakest for populations of states with core penetrating 
valence electron orbitals and we expect spin system breakdown for high nl states progressively 
at lower n and I for increasing ion charge z- Thus the above assumptions are appropriate 
for light element ions with a more elaborate bundle-{J p)nl j model more suited to heavy 
element ions. The latter will be the subject of a separate paper. In the hnndle-nS model, 
only equilibrium populations of complete n-shells for a given parent and spin system need be 
evaluated, which are the solutions of the statistical balance equations 



n-1 



n"=no 

CO 

-{T^[Neq':l,,+N,qit„,,] (32) 

n'=n+l 
n-\ 

+ J][A„^n"+N4:l„„+N,ql^l„„] 

«"=no 

V' 

is the population of the parent ion Ay~' , is the free electron density and A^^, the 
free proton density, A is the usual Einstein coefficient, q^"^ and q'^''^ denote coUisional rate 
coefficients due to electrons and protons, alP, a[f' and aj,^' denote radiative, dielectronic and 
three-body recombination coefficients respectively, A"^^^^^, denotes secondary autoionisation 

and q^nl^g and q\^}^g denote collisional ionisation rate coefficients due to electrons and protons 
respectively. cos,,s is a spin weight factor (see equationl^below). For complex ions, there are 
separate systems of equations for each parent and for up to two spin systems (in L-S coupling) 
built on each parent and one such equation for each value of n from the lowest allowed n-shell 
«o for the parent/spin system to oo. The number of spin systems is labelled A^.,,., and the lowest 

if) 

allowed n-shell no = n,) , the lowest accessible shell by recombination except for doublets 
built on the He-like \s2s parent. In general bare nuclei of other elements are effective ion 
projectiles along with protons. We use the word protons here to represent mean Ze// ions with 
suitably scaled collisional rate coefficients. These equations are analogous to the equations for 
hydrogen and coupling-independent expressions may be used for the main « — > «' coefficients 
providing a suitable spin system weight factor 

= 2(25771) ^^^^ 
is introduced. Table [2 summarises these various parameters for first period iso-electronic 
sequences up to fluorine/neon. There are a number of issues. 
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2.1.1. b-factors and lowest levels For hydrogenic ions it was advantageous to write the 
statistical equations in terms of Saha-Boltzmann deviation factors. This remains true for 
complex ions but the definition must be generalised. The deviation by „s is defined by 



2t \3/2 



exp(/v.„s /kTe)bv,„s ■ (34) 



kTe ) 2u)y 

That is by„s is specified with respect to the parent ion state A-^~\ and the ionisation potential 
/y„s is also referred to that parent, oq is the Bohr radius. Note that Wv,„s - LJs,,,sfi^^v 
where lOy — (25 v + l)(2Lv + 1) is the parent statistical weight. It is convenient to introduce 
Cv,ns - bv,ns - 1 ^nd scaled temperatures and densities 
'kTAl_ 

5 l^TialN, 

with similar forms for the proton temperature and density, a is the fine structure constant. In 
these terms the statistical balance equations become particularly suitable for calculation and 
resulting 963? coefficients (z-scaled), on a Oelpe grid, can sustain interpolation of moderate 
precision within an iso-electronic sequence. 

The parent/spin system model does not distinguish recombined metastables within the 
same spin system. For example, consider the 2s^2p parent in a B-like like ion recombining 
into the singlet system of the C-like ion. The parent/spin system has two metastables, 
2s^2p^ 'D and 2s^2p^ '5 . We assign the effective ground state as the lowest energy metastable 
(2s^2p^ 'D state in this example) and assign statistically weighted n-shell sub-populations to 
resolve between the recombined metastables. There is often a substantial difference between 
the quantum defect of the lowest state itself and the mean quantum defect arising from a 
weighted average of the lowest n-shell terms. These choices improve the stepwise part of 
the collisional-radiative ionisation coefficient within the bundle-nS' picture at low density and 
allow the bundle-nS' model to stand alone, albeit at some reduction in accuracy. A precise 
solution, as adopted in this paper, is given by the condensation/projection/expansion matrix 
transfer of the hnndle-nS model onto the low level solution which distinguishes explicitly 
between the LS terms in the lowest few quantum shells from the beginning (see section IT^ . 



2.7.2. Auto-ionisation and alternate parents Let the metastable state which represents 
the lowest level of recombined parent/spin system v, nS be labelled by p. With parent v, 
from a stepwise collisional-radiative point of view it is convenient to refer to v, nS as an 
intermediate state system and label it (v) '^■^"^'^n. Thus a recombination/stepwise/cascade 
pathway may progress from a parent v through the intermediate state y, nS ending on the 
metastable p with the effective SC3? recombination coefficient written as yiC'D,,^(v) (2s+n„;p and 
likewise a stepwise excitation/ionisation pathway may progress from lowest metastable p via 
intermediate state v, nS to final ionised ion metastable v with the effective SC3? ionisation 
coefficient §C'Dp^(v) (2s+i)„;v. The population calculation for a given pathway involves the 
excited state population structure connecting the recombining parent v and metastable p. 
However, an alternative parent, which is not the intermediate state parent, can be populated 
by autoionisation of excited states above the auto-ionisation threshold. That is, if the parent v 
is a metastable and so there exists a lower lying ground state (or possibly other metastable) of 
the parent, say v' . The excited state populations of the v,nS system must include such auto- 
ionisation processes. Above the auto-ionisation threshold, and at low electron densities, auto- 
ionisation is the dominant loss mechanism. However, the bundle-n auto-ionisation transition 
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Figure 3. Bundle-«5 population structure for O^^ recombining from and ionising to 0^^ 
at Tg = 1.7 X K. (a) hf^^i is) 'n factors for doublet bundle-n populations built on the 
ground parent 2s- '5 term. The enhanced b-factors at n ~ 10 - 300 is due to dielectronic 
recombination. Curves show the progressive suppression of the high populations as the 
electron density increases, (b) b^2„p ip) 2„ factors for doublet bundle-n populations built on 
the metastable parent 2.v2p term. The abrupt transition to underpopulations at « ~ 4 is due 
to secondary Auger transition {LS coupling allowed breakdown type) to the ' 5 parent system. 
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probabilities scale as rT^ and are independent of electron density, whereas direct ionisation 
loss rates scale as n"* and vary directly with electron density. For high n-shells, direct 
ionisation is the dominant loss process. As the electron density increases, direct ionisation 
becomes the dominant loss process for all n-shells and auto-ionisation is quenched. The 
inclusion of auto-ionisation transition probabilities in the statistical balance equations leads 
to dramatic changes in the population structure as shown in figures [3^, b. In the example, 
the populations of C^^{2s^ '5) are built upon a ground state parent so no auto-ionisation 
pathways are accessible from the excited states. If the populations are expressed in terms of 
the Saha-Boltzmann b-factors, then they show strong overpopulation of the high n-shells at 
low electron density due to dielectronic recombination. As the electron density is increased, 
the dielectronic recombination influence becomes less due to ionisation of the high n-shell 
populations with all the b-factors tending to 1 . This behaviour is typical for a recombined 
system built on a ground state parent. By contrast, the populations of C^^(2s2p ^P) ^n, built 
upon an excited metastable parent, show powerful depopulation above the auto-ionisation 
threshold. As the electron density increases, the direct ionising collisions from the excited 
states compete more strongly and the b-factors tend to 1 . 

The definitions of generalised collisional-radiative ionisation and recombination 
coefficients are still relevant but give the recombination rate from v to p and the total loss rate 
from p (i.e. with no resolution of final parent after ionisation) respectively. The recombination 
coefficient is already correctly parent/metastable resolved and needs no further adjustment. 
The correct parent resolved (partial) ionisation coefficients §C'Dp^^y^ <^-s*^>ny derived by 
constructing the loss vector from each level. Because the model considers the excited state 
populations Nv,„s , direct ionisation only populates Ay^' . The alternative parents are populated 
by auto-ionisation. The parent resolved loss vectors are thus given by 



lonisation pathways and the expected (partial) 9631 ionisation coefficients for C^' -H e 
C^^ + e + e are summarised in table |2] From the point of view of recombination from v via 
intermediate states v, nS towards p, autoionisation allows exit into the alternative parent v' 
before p is reached. This gives the new parent cross coupling coefficient XCDy^^y-, (2s+i)„ y,. 
Note that in applications, it is the coefficients summed over intermediate states which are 
required as 



Behaviours are illustrated in figures 0^, b. 
2.2. The low-level + projection model 

Consider now the set of low levels of an ion. For 9 63? studies, we are concerned with complete 
sets of low levels associated with the valence electron in the ground and excited n-shells. The 
span of the low levels is Amqi - n\ - no with no denoting the ground complex valence n-shell 
and «i the highest resolved n-shell. For light elements and spectroscopy which extends up to 
the visible range, we seek Anoi = 1 at minimum and 2 preferably. 

The solution for the populations and effective coefficients follow the theory of section 
11.2. H and is numerically straightforward, carried out in the population representation (called 
the /^-representation), rather than ^-factor or c-factor representations (see equation 1?^ and 
the following lines). The use of the c-factor representation, where c - b - 1, is essential for 




(36) 




(37) 



v,S 
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Figure 4. 



Bundle-«S generalised collisional-radiative coefficients. 

-+2 



Tg is the electron 

temperature and Ne the electron density, (a) Cp' + e — > Cj" + e + e for difl'erent initial 
metastable p, final states v and intermediate parent/spin systems. Curves are labelled as 
-> (v') v) where p = \= (2s^2p ^P),p = 2 = (ls2p^ v,v' = 1 = (2*^ 's) 

and v,v' = 2 = {2s2p ^P). (b) 0*^[2s2p^ '•P} + e ^ 0+^[2s^2p ^P] + e parent cross- 
coupling coefficients. Both 0^~[(2.s2p^ ^P) and 0*\(2s2p^ "^P) ^n] recombined systems 
are included. The autoionisation pathway for the former is allowed in LS-coupling, but for 
the latter proceeds only in intermediate-coupling (see section IT^ . Such coefficients follow 
broadly the behaviour of recombination coefficients. 
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Final parent Intermediate Initial SC3?ionis. Type 
metastable system metast. coefft. 



1 


(2s^ 'Syn 


1 


862)1^(1) 2„j 


d 


1 


{2s2p ^P) 


1 


SCDi_(2) 2„j 


e-a 


2 




1 


§C23i^(2) 2„_2 


d 


1 


(2s2p V) 4n 


2 


SCD2-»(2) "n.l 


e-a (IC) 


2 




2 


SCD2^(2) in,2 


d 



Table 2. Excited state structures and partial ionisation coefficients calculated for the Be- 
B series. lonisation coefficient nomenclature is SCCD^,_,(j,j {25+i)„ yf where p indexes the lowest 
metastable of the parent/spin system, with y the parent of the spin system and v' indexes the 
final metastable state, 'd' denotes direct lonisation and 'e-a' denotes inner shell excitation- 
autoionisation contributions. 



cancellation error avoidance in computation of the very high level population structure, but is 
not necessary for the low levels. The construction of the Ctj and r,_v coefficients is by spline 
interpolation in electron temperature of data extracted from archives. Expansion of indirect 
projection data and its amalgamation with the resolved direct data is by weight matrices as: 



■:v,S 



C_ /^dir , \ ^ , c /^indir 
v,S 

r r dir , \ ^ , X r indir 

Li - L. + 2j W,:v,s^,:v,5 



where / e n and j e n' . Introducing the term statistical weight fractions as 
_ (25 ,■ + 1)(2L,. -H 1) 

"2:,e„(25, + l)(2L, + l) ^'^^ 

then 

^/;:y,S = '^S,.S^^/:«:v,S for U + u' 





s = ws, 


.,s 




K.v,S 




r.v.S 


for n = n' 











and cJSy.s is as defined in equation|33l Table|3]illustrates the C,y weighting. 

In the deduction of spectral emission coefficients between low levels there can be some 
confusion. Following the definition of equation |2] in the resolved low level picture, the 
emission coefficient is referred to a particular metastable. If metastables are neglected, so 
that there is only a ground state and all other levels are viewed as excited, the reference is 
to the ground state. On the other hand, a metastable treated as an ordinary excited level will 
have a quasistatic population comparable to that of the ground so that 2, M = N,o, + N\ and 
J'£C/o/j'^A: = y^^ij^kNxINioi and this is the coefficient which should be used with a 'stage- 
to-stage' (that is un-generalised picture) lonisation balance. Figure|5lillustrates the behaviour 
of low level populations. The graph is of the parameter 3^*"^' from eauationl20l 
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'ndex 


Term 


Spin 


Parent 


Shell 


n-shell weights 
n=2 n=3 


1 


2s^2p "-P 


2 


1 


2 


0.25 


1 


2s^2p 


2 


2 


2 


0.25 


2 


2s2p^ '^P 


4 


2 


2 


0.786 - 


3 


2s2p^ 


2 


1 


2 


0.45 


3 


2s2p^ 


2 


2 


2 


0.45 


4 


2s2p^ 


2 


1 


2 


0.05 


4 


2s2p^ ^-S 


2 


2 


2 


0.05 


5 


2s2p^ 2p 


2 


2 


2 


0.25 


5 


2s2p^ 2p 


2 


2 


2 


0.25 


6 


2^23525 


2 


1 


3 


0.0667 


6 


2^23^25 


2 


2 


3 


0.0667 


7 


2^23^ 2p 


2 


1 


3 


0.3333 


7 


2^23^ 2p 


2 


2 


3 


0.3333 


8 


2p3 4p 


4 


2 


3 


0.214 - 


9 


2s^3d^D 


2 


1 


3 


0.6 


9 


2s^3d 2d 


2 


2 


3 


0.6 


10 


2s2p3s^P 


4 


2 


3 


0.1089 


11 


2s2p3p 


4 


2 


3 


0.1881 


12 


2s2p3p '^S 


4 


2 


3 


0.0297 


13 


2s2p3p ^P 


4 


2 


3 


0.1089 


14 


2s2p3d^F 


4 


2 


3 


0.2673 


15 


2s2p3d^D 


4 


2 


3 


0.1881 


16 


2s2p3d'^P 


4 


2 


3 


0.1089 



Table 3. Fractionising of bundle-n5 matrix projection onto low levels for B-lilce ions witli 
two resolved n-sliells. 



2.3. Specific reactions 

Energy levels for high bundle-n5 levels and their A-values, Maxwell averaged collision 
strengths, radiative recombination coefficients, dielectronic recombination coefficients and 
ionisation coefficients are generated from a range of parametric formulae and approximations 
described in earlier works (Burgess and Summers, 1976; Summers, 1977; Summers and 
Hooper, 1983; Burgess and Summers, 1987; Summers, 2004). High quality specific data 
when available are acccessed from archives and substituted for the default values. This is a 
systematic procedure for dielectronic coefficients (see sections |2 . 3 . 1 l and Islbelow) . 

For low levels, a complete basis of intermediate coupling energy level, A-value and 
scaled Born approximation collision data, spanning the principal quantum shell range no ^ 
n < «! is generated automatically using the Cowan (1981) or Autostructure (Badnell, 
1986) procedures. This is called our baseline calculation. These data are merged with 
more restricted (in level coverage) but similarly organised highest quality data from archives 
where available (e.g. R-matrix data such as Ramsbottom et al (1995)). The data collection 
is compressed by appropriate summing and averaging to form a complete LS term basis 
and augmented with comprehensive high quality LS resolved dielectronic recombination, 
radiative recombination and collisional ionisation coefficients mapped from archives (see 
sec tions [2 . 3 . 1 1 12 . 3^ and l5l below) . The radiative data has its origin in the work of Burgess 
and Summers (1987). The final data collection is called a specific ion file (ADAS data format 
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Figure 5. Low level population structure of Ne^^. The excitation part driven by the \s^2s 
ground state is denoted by Ni . Terms are labelled as 2: 2p ^P, 3: 3s ^S, 4: 3p ^P, 5: 3d ~D, 
6: 4.V , 7: 4p ^P, 8: Ad 'D. At low density, the population ratios Ni/NeNi become flat 
tending to the coronal values. At high density the population ratios decrease inversely with 
the electron density, Ng, tending to the Saha-Boltzmann values. The intermediate region is the 
collisional-radiative regime. 



adf04). The detailed content is examined in section|3 Details of state selective dielectronic 
recombination and ionisation coefficients are given in the following sub-sections. 

2.3.1. State-selective dielectronic recombination State selective dielectronic recombination 
coefficients are required to all resolved low levels and to all bundle-n5 shells for the various 
initial and intermediate state metastable parents v for modelling. These are very 

extensive data and have been prepared for the present 363? work through an associated 
international 'DR Project' summarised in Badnell et al (2003) [hereafter called DR - paper I\ 
and detailed in subsequent papers of the series . Based on the independent particle, isolated 
resonance, distorted wave (IPIRDW) approximation, the partial dielectronic recombination 
rate coefficient a^''^,. from an initial metastable state v into a resolved final state / of an ion A^^- 
is given by 



(d) 

a . 



A 2, \3/2 



(41) 



where ojpj is the statistical weight of the (A^+ l)-electron doubly-excited resonance state j, (±>y 
is the statistical weight of the A^-electron target state and the autoionization {A^) and radiative 
(A') rates are in inverse seconds. The suffix p is used here to denote a parent ion state. is the 
energy of the continuum electron (with angular momentum /), which is fixed by the position 



H.P. Summers et al. lonisation state. I. 



20 



of the resonances. Note that the parent states p are excited, that is they exclude the metastable 
parents v'. The code autostructure (Badnell, 1986; Badnell and Pindzola, 1989; Badnell, 
1997) was used to calculate multi-configuration LS and intermediate coupling energy levels 
and rates within the IPIRDW approximation. The code makes use both of non-relativistic and 
semi-relativistic wavefunctions (Pindzola and Badnell, 1990) and is efficient and accurate for 
both the resolved low level and high-n shell problems. Lookup tables (see section|5l adf09) 
are prepared comprising state selective recombination coefficients at a standard set of z-scaled 
temperatures, for each metastable parent, to all LS resolved terms of the recombined ion with 
valence electron up to n-shell n*''' > n\ (normally n"^^ - 5) and to bundle-n5 levels of a 
representative set of n-shells (usually spanning no to 999). These bundle-n5 coefficients are 
simple sums over orbital states / e n and so apply at zero density. This provides an extensive, 
but still economical, tabulation. 




Figure 6. O dielectronic recombination, bp ,,; factors as a function of outer electron orbital 
angular momentum / for doubly-excited states of O^^ relative to 0*+ 2s' ' S for p = 2s2p ' P 
and 2s3p 'P, n = 20, T, = 10^ K, = Np, Z,ff = 1. Cases: 1. = lO'" cm^^; 2. = lO'^ 
cm"^^; 3. Nc = 10'^ cm"^; 4. Nc = lO''' cm"'; 5. = lO'^ cm"'. Note the alternative Auger 
channel reduction for the p = 2s3p 'P graphs. 



In DR paper-I, we introduced an associated code, the 'Burgess-Bethe general program' 
BBGP BBGP is used here as a support function in a model for the /-redistribution of doubly- 
excited states which provide a correction to the accurate, but unredistributed, dielectronic data. 
The redistributed data (regenerated in the same adf09 format) are normalised to IPIRDW at 
zero density. The procedure is similar to that for singly excited systems. For LS' -averaged 
levels, the number densities expressed in terms of their deviations, bp^„i from Saha-Boltzmann, 
and referred to the initial parent v', are given by 



3/2 



p,nl • 



(42) 



Thus, in the BBGP zero-density limit, with only resonant capture from the v parent balanced 
by Auger breakup and radiative stabilization back to the same parent, we have 

y, 4^ 

^1 nl^v'lt'V 
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p,nl 



In the extended BBGP model, we include resonant-capture from initial metastables other 
than the ground, dipole-allowed collisional redistribution between adjacent doubly-excited 
/-substates of the same n by secondary ion- and electron-impact, and losses by 'alternate' 
Auger break-up and parent radiative transition pathways. The population equations for the 
/-substates of a doubly-excited n-shell become 

V/'=/±i /'=;±i 
p-i i+i p-i \ 

p, = ll'=l-l p,=l 

M l+l P 
v=l /'=/-! Pl=P+l 

q'^ denotes resonance capture coefficients. M denotes the number of parent metastables 
which are the starting point for resonance capture whereas P denotes the number of true 
excited parent states. Ion impact redistributive collisions are effective and are represented 
in the equations as Zeff ion contributions. The density corrected bundle-«5 recombination 
coefficients are then 

(dJPIRDW),^ . WPIRDW),^ _U_^^^y^^nl_^_e>_ 

Details are given in DR paper-I. Figure |6l shows the effects of redistribution of the doubly 
excited bp^„i factors in the recombination of O^'*. 

Autoionisation rate calculations in LS coupling exclude breakup of non-ground 
metastable parent based 'singly excited' systems by spin change. Such breakup only occurs 
in intermediate coupling but cannot be ignored for a correct population assessment even in 
light element ions. These rates are computed in separate intermediate coupling dielectronic 
calculations. These rates are main contributors to the parent metastable cross-coupling 
identified earlier (cf. the 0^^[{2s2p^ '^P) ^n] intermediate states in figure leading back 
to the 0^^(2s^2p ^P) parent). 



2.3.2. State-selective ionisation Direct ionisation coefficients for excited n-shells in bundle- 
nS population modelling are evaluated in the exchange-classical-impact parameter (EClP) 
approximation (Burgess and Summers, 1976). The method which merges a symmetrised 
classical binary encounter model with an impact parameter model for distant encounters is 
relatively simple and of moderate precision, but has a demonstrated consistency and reliability 
for excited n-shells in comparison with more eleborate methods. There is a higher precision 
requirement for state selective ionisation from metastable and low levels. Such ionisation 
includes direct and excitation/autoionisation parts, but the latter only through true excited 
parents. Stepwise ionisation is handled in the coUisional-radiative population models. There 
are extensive ionisation cross-section measurements, but these are in general unresolved and 
so are of value principally for renormalising of resolved theoretical methods. The most 
powerful theoretical methods used for excitation, namely R-matrix with pseudostates (RMPS) 
(Bartschat and Bray, 1996; Ballance et al , 2003) and convergent-close-coupling (CCC) (Bray 
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and Stelbovics, 1993), have some capability for ionisation, but at this stage are limited to 
very few electron systems as is the time-dependent-close-coupling (TDCC) (Pindzola and 
Robicheaux, 1996; Colgan et al , 2003) ionisation method. For SC3^ calculations, we have 
relied on the procedures of Summers and Hooper (1983) for initial and final state resolution 
of total ionisation rate coefficients and on the distorted wave approximation. The distorted 
wave method is our main method for extended SC3^ studies for many elements and most 
stages of ionisation. RMPS and TDCC studies are directed mostly at the neutral and near 
neutral systems where the distorted wave method is least reliable (Loch et al , 2005). We 
use the configuration average distorted iCADW) wave approach of Pindzola, Griflin and 
Bottcher (1986). It has reasonable economy of computation while allowing access to complex, 
multi-electron ions, highly excited states, excitation/autoionisation and radiation damping. It 
expresses the configuration averaged excitation cross-section for 

(n,hr^\n2hr-%li ^ (nxhfHmhf'lflf (46) 

as 

87r 

O-excit = + l)(4/2 + 3 - g'2) 

k^kf 

Yjik + mif + mif-M (47) 

and the configuration averaged ionisation cross-section for 

(nihr^%li ^ (mhf'kMflf (48) 



as 



O-ionis = ■- (qi + 1) 



(2Z,- + 1)(2/, + 1)(2Z/ + l)M,f,u (49) 

liJeJf 

where M2/;i, are squared two-body Coulomb matrix elements, o" denotes the average cross- 
section and kj,kf and k^ denote average initial, final and ejected electron momenta respectively 
in the configuration average picture. The configuration average direct ionisation cross- 
sections and rate coefficients may be unbundled back to resolved form using angular factors 
obtained by Sampson and Zhang (1988). Note that the excitations described here are to auto- 
ionising levels and so resolved Auger yields may be used which are the same as those in the 
dielectronic calculations of sec tion 17. 3 . 1 1 above . In fact the ionisation coefficient calculation 
results are structured and archived in ADAS (format adf23) in a manner very similar to state 
selective dielectronic recombination. Extensive studies have been carried out on light and 
heavy elements (Colgan et al , 2003; Loch et al , 2003). 



3. Fundamental atomic data for low levels of ions 



For bundle-n5 modelling, the expected fundamental data precision is < 30% for excitation and 
ionising collisional rate coefficients, < 10% for A-values and state selective recombination 
coefficients and < 1% for energies. 

For low level modelling many sources are used. A rating is given for the classes 
of fundamental data for the ions of carbon, oxygen and neon in table |4] which is based 
on the following considerations. For energy levels, categories are a spectroscopic, b < 
0.5% and c < 1.0%. Category c is anticipated from ab initio multi-configuration structure 
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Table 4. Fundamental data precision categorising for the ions of carbon, oxygen and neon. 
The definition of the categories and justification of the categories for Ne** are given in the 
text. 



calculations, b from such calculations with extended optimising and a reflects direct inclusion 
of experimental energies from reference sources. For A-values, categories are a < 5%, b 

< 10% and c < 25%. Category c is anticipated from our baseline calculations, b from 
optimised multi-configuration structure calculations with extended optimising and a from 
specific studies in the literature. For electron impact Maxwell averaged collision strengths, 
T, the categories are a < 10%, b < 20% and c < 35%. Category c is from our baseline 
calculations, b from distorted wave calculations and a from specific R-matrix calculations, 
equivalent methods or experiment. For radiative recombination, the categories are a < 5%, b 

< 10% and c < 50%. Category c is from scaled methods using hydrogenic matrix elements, b 
from distorted wave one-electron wave functions in an optimised potential using spectroscopic 
energies and a from specific R-matrix calculations and experiment. Category b is the baseline 
in ADAS. For dielectronic recombination, the categories are a < 20%, b < 30% and c < 
45%. Category c is from BBGP approximations, b from LS-coupled calculations using 
Autostructure from the DR Project and a from IC-coupled calculations using Autostructure 
with parent and lowest resonance energy level adjustments from the DR Project. Category b is 
the baseline in ADAS. It is to be noted that the theoretical relative precision which is consistent 
with the variation between the three categories is 15% better but dielectronic recombination 
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comparisons with experiment still show unexplained discrepancies at the 20% level so the 
present categories are safe. For ionisation, the categories are a < 10%, b < 25% and c < 40%. 
Category c is from ECIP approximations, b from configuration average with angular factor 
term resolution and a from RMPS, TDCC calculations and experiment. Category b is the 
baseline in ADAS. 

As indicated in section 12.31 these various data are assembled in an adf04 file which is 
sufficient to support the primary low-level population calculation. The tabulations are at a set 
of temperatures arising from a fixed set of z-scaled temperatures (see section I2.1.1> which 
spans the full range to asymptotic regions of reaction data. Collision data are converted 
to these ranges using C-plots (Burgess and Tully (1992) and this procedure also flags data 
errors or queries. The precision of the specific ion file determines the achievable precision 
of all derived populations, emissivities and collisional-radiative coefficients (Whiteford et al , 
2005; O'Mullane et al , 2005). The assembling of data in the adf04 is systematic and orderly. 
A comment section at the end of the file details the assembly steps, implementer, codes 
used and dates. This includes baseline and supplementation files, merging, LS compression, 
dielectronic recombination data inclusion etc. Also there is extended detail of orginal data 
sources and a history of updates. A given adf04 file represents a 'snapshot' of the state 
of available knowledge at the time. It is subject to periodic review and ADAS codes (see 
section|5} are designed to enable easy reprocessing of all derived data following fundamental 
data update. The grading for each ion given in table 0] is justified and supported by the 
comments from its adf04 file. A summary from Ne^^ is given in the following paragraphs 
in illustration. For Ne^^, the low levels span 44 terms, including up to the n - 5 shell 
built on the parent \s^2s '5 and n = 3 shell built on the \s^2p '5. Only the \s^2s '5 
parent is treated as a metastable from the SC3^ point of view. The intermediate coupled 
baseline data set is copmm#10-ic#ne6.dat with preferred supplementary energy A-value and 
T data merged from copjl#be-ic#ne6.dat. lonisation potentials and energy levels are from the 
NIST Standard Reference Database apart from 2p3p '5 (Kelly, 1987) and 2p3d levels by 
(Ramsbottom et al , 1995). The categorisation is a. A-values were drawn from the Opacity 
Project( Opacity Project Team, 1996; Tully et al , 1991 ) as justified (for N IV and O V) by 
Wiese, Fuhr and Deters (1996) and supported/adjusted by < 3% by Fleming et al (1996a,b), 
Jonsson et al (1998), Froese Fischer, Godefeid and Olsen (1997), Froese Fischer, Gaigalas 
and Godefried(1997), Nussbaumer and Storey (1979), Sampson, Goett and Clark (1984) and 
Ramsbottom et al (1995). Category b is safe. Ys are taken from Ramsbottom et al (1995). 
These are for a 26 LS eigenstate multichannel R-matrix calculation. The categories assigned 
are a and b. Radiative and dielectronic recombination and ionisation all follow the baselines 
in ADAS, that is categories b although the summed and averaged ionisation rate coefficient 
(over metastables) is normalised to experiment. 

4. Illustrative results 

Figure0shows IP£C coefficients for the C II 858 A spectrum line. The coefficients depend on 
both electron temperature and electron density in general. A common practice in spectral 
analysis is to observe principally the strongest resonance line of an ion. Such emission 
is driven largely from the ground state (figure Et) and because of the large A-value, the 
density sensitivity occurs at relatively high density. Thus such resonance line emission at 
moderate to low densities mostly reflects temperature and the distribution of the ionisation 
stages. Comparison in near equilibrium ionisation balance plasmas of line ratios from the 
same ion, by contrast, is mostly directed at electron density and relies on the presence of 
metastables and spin changing coUisional processes to confer the sensitivity. As shown 
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5.0"^ F (a) 




Figure 7. Excitation and recombination plioton emissivity functions, CP£Cs, driven by tlie 
ground metastables, are shown for the C 11 858 A spectrum line. The coefficients depend 
on both electron temperature and electron density. The (exc) pail decreases at high electron 
density due principally to stepwise lonisation losses from the upper level of the transition. 
The behaviour of the (red) part shows both suppression of dielectronic recombination at 
moderate density and then enhancement due to three-body recombination at high density and 
low temperature. 
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Figure 8. Generalised collisional-radiative recombination coefficients for + e — » 0^ (a) 

AeD(2s^ 'S -» 2s22p -P) (b) Ae'D(2i'^ '5 -> 2s2p2 "p) (c) yieXl(2s2p ^ 2s2p2 *P) 
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earlier, the balance of the dominant ground and metastable populations is disturbed in dynamic 
plasmas and so density sensitivity may be modified by the dynamic state. The distinction of 
the metastable driven TfiC in SC3^ modelling, as illustrated in figure03, allows more complete 
diagnostic study and the possibility of separation of the two effects. In strongly recombining 
plasmas (most commonly photo-ionised astrophysical plasmas) the direct contribution of 
recombination to the emission may dominate the excitation part. The ^PfiC*"^'*, as illustrated 
in figure 0;, is then required. In the fusion context, multi-chordal spectral observations are 
important for the study of impurity transport especially near sources. Visible and quartz UV 
observations are convenient and this places a requirement for TfiCs from higher quantum 
shells. It is this primarily which defines the span of our low-levels for population modelling. 
Emission from higher n-shells is significantly aff'ected by cascading from yet higher levels. 
The full machinery of projection as described in section 12.11 is necessary for our global 
ambition of 20% precision for emissivity coefficients. 

Figure illustrates the 9631 recombination coefficients. At low electron density, 
radiative and dielectronic recombination dominate. For capture from metastables, alternate 
Auger branching can largely suppress the dielectronic part of the surfaces. Thus figure [8^ 
shows the characteristic exponential rise at the temperature for excitation of the main parent 

—3/2 

(core) transition of dielectronic recombination and then the subsequent fall-off' (as ~ ), 
in contrast with |8j;. At moderate densities, suppression of the high «-shell populations, 
principally populated by dielectronic recombination, through re-ionisation occurs and the 
coefficient falls in the dielectronic recombination region. At very high density, three-body 
recombination becomes effective, preferentially beginning at the lower electron temperatures. 
It is evident that recombination is less effective from the metastable at relevant ionisation 
balance temperatures. Models which ignore the role of capture from the metastable parent 
(which may be the dominant population) can lead to substantial errors in recombination 
coefficients, while simply excluding all capture from the metastable cannot deliver the 
precision sought for current modelling. For light elements in astrophysical plasmas, the zero- 
density coronal assumption for the recombination coefficient is still frequently made. This 
cannot be justified even at solar coronal densities. 

Figure |5] illustrates the QCJi ionisation coefficients. At low electron density, the 
coefficient is dominated by direct ionisation, including excitation/auto-ionisation, from the 
driving metastable. Relatively high electron densities are required before the stepwise 
contribution begins. It is primarily the excitation to, and then further excitation and ionisation 
from, the first excited levels which controls this. The ground and metastable resolved 
coefficients both show the same broad behaviour as the usual (stage to stage) collisonal- 
radiative coefficient, tending to a finite limit at very high density. It is to be noted that the 
XCD coefficients are required to be able to construct a meaningful stage to stage collisional- 
radiative ionisation coefficient from the generalised progenitors. It remains the case that most 
plasma modelling (certainly in the fusion area) is not adjusted to the use of the generalised 
coefficients as source terms. Reconstruction of stage to stage source terms (at the price of a 
reduction in modelling accuracy) from the generalised coefficients is still is a requirement and 
is addressed more fully in section|5] 

The generalised coefficients may be used to establish the equilibrium ionisation balance 
for an element in which the dominant ground and metastable populations are distinguished, 
that is the fractional abundances 




: cr = l,---,M,;z = 0, •■■,zo 



(50) 



Nltot] ^ 
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where M; is the number of metastables for ionisation stage z and 

-'0 Zo 

^[M^J^W^J^^W (51) 

z=0 2=0 o-=0 

in equilibrium. Writing N}~^ for the vector of populations n]^^, the equilibrium population 
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GCR (zero density) gcr 
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1 10 100 1000 
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Figure 10. Equilibrium ionisation balance of oxygen (a) Equilibrium fractional abundances 
of metastables from the generalised collisonal-radiative (b) The total radiated power function 
showing the components driven by the different metastables (c) The Qiyi function for the O II 
629.7 A line showing the contributions parts from the different metastables of the z-times and 
(z + l)-times ionised ions. 
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fractions are obtained from solution of the matrix equations. 



g[0,0] 






g[l,l] 





0] 





g[2,2] 



1] 



^[2] 



= 0. 



(52) 



equil 



These in turn may be combined with the CPZLT and CPD?® to obtain the equilibrium radiated 
power loss function for the element as 



p{tot] ^ ^ p{z. 

z=0 



equil 



2=0 o-=0 



(53) 



equil 



The equilibrium fractional abundances and equilibrium radiated power function are illustrated 
in figure^! It is useful at this point to draw attention to emission functions which combine 
emission coefficients with equilibrium fractional abundances. They are commonly used in 
differential emission measure analysis of the solar atmosphere (Lanzafame et al , 2002) 
where they are called G(Te) functions. In solar astrophysics, it is assumed that the G{Te) 
are functions of the single parameter (either from a zero-density coronal approximation or 
by specification at fixed density or pressure) and usually the abundance of hydrogen relative to 
electrons Nh/N,, is incorporated in the definition. For finite density plasmas, in the generalised 
collisonal-radiative picture, we define QTN functions, parameterised by and A^^, as 



■'equil 



(54) 



/ equil 



These are strongly peaked functions in T^. The precision of the present modelling and data, 
including the full density dependence, is in principle sufficient to allow bivariate differential 
emission measure analysis (Judge et al , 1997) although this remains to be carried out. The 
general behaviour of a STN function is illustrated in figure [Tot. 



5. Computations and archiving derived data for applications 

The organisation of the main calculations and data flow is shown in figure^] The calculations 
executed for the paper have been implemented in general purpose codes and attention has 
been given to the precise specification of all data sets and the machinery for accessing and 
manipulating them. This includes initial data, intermediate and driver data as well as the final 
products and follows ADAS Project practice. Population structure modelling requires 
an initial input dataset of energy levels, transition probabilities and collisional rates of ADAS 
data format adf04 which is complete for an appropriate designated set of low levels. For 
the light elements, best available data were assembled and verified as described in section|3] 
These data are substitutes for more moderate quality, but complete, baseline data prepared 
automatically (see GCR - paper III). The adf04 files are required for LS-terms. In practice, 
we find it most suitable to prepare the data for LSJ levels and then bundle back to terms. For 
SC3^ modelling, state selective recombination and ionisation coefficients must be added to 
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Figure 11. Schematic of computational steps in production of generalised collisonal-radiative 
data. Data sets are shown by circular outline and program elements by rectangles. The key 
members have been incorporated and assigned names in the ADAS system. 



form the fully specified adf04 file. The LS-coupled dielectronic coefficients are mapped in 
from very large comprehensive tabulations of data format adjD9 prepared as part of the DR 
Project (Badnell et al ,2003). The resulting adf04 files are the ADAS preferred data sets and 
are available for elements helium to neon. 

The two population codes, called ADAS204 and ADAS208, work together. ADAS204 
is the bundle-n5 model. Data on the metastable parent structure, quantum defects, auto- 
ionisation thresholds and autoionsation rates are required. These data may be extracted from 
the adf04 files and it has been helpful to make the driver dataset preparation automatic. High 
quality shell selective dielectronic data are essential and this is part of the provision in the 
adf09 files described in the previous paragraph. ADAS204 provides complete population 
solutions, but extracts from these solutions the condensed influence on the low n-shells, as 
projection matrices, for connection with the calculations of ADAS208. It is to be noted that 
the main on-going development is refinement of atomic collision rates between the key low 
levels. The projection matrices are not subject to frequent change and so are suitable for long- 
term archiving (adfl7). ADAS208 is the low-level resolved population model which delivers 
the final data for application. It draws its key data from the fully configured adf04 file and 
supplements these with projection data. The evaluation of the population structure takes place 
at an extended set of z-scaled electron temperatures and densities (see section O and this 
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means that the resulting coefficient data are suited to interpolation along iso-electronic 
sequences. Thus the initial tabulation of SCIR coefficients is in iso-electronic datasets. It is 
convenient to implement the gathering and mapping from iso-electronic to iso-nuclear in a 
separate step, which also supports the merging back to the unresolved stage-to- stage picture 
if required. The production of CPfiC and coefficients is directly to iso-nuclear oriented 
collections. The number of JSCs from the population calculations can, in principle, be very 
large. We restrict these by a threshold magnitude and to particular important spectral regions. 
It is straightforward to rerun ADAS208 to generate CPS 6 or SXCB coefficients alone in spectral 
intervals of one's choice. The separation of the ADAS204 and ADAS208 tasks and the ease 
of modifying data within an adjD4 file means that 'what-if ' studies on the sensitivity of the 
derived data to fundamental data uncertainly can readily carried out. Special ADAS codes 
enable detailed study of cumulative error and dominating sources of uncertainty such that 
error surfaces, for example for a TEC as a function of electron temperature and density, may 
be generated. Such error (uncertainty) analysis for theoretical derived coefficients and its 
utilisation in the confrontation with diagnostic experiments is the subject of a separate work 
(see O'MuUane etal , 2005). 

6. Conclusions 

The requirements for precise modelling of spectral emission and the relating of ionisation 
stages in thermal plasmas have been considered. ColUsional-radiative methodologies have 
been developed and extended to enable the fuU role of metastables to be realised, so that 
this generalised (SC^K) picture applies to most dynamically evolving plasmas occurring in 
magnetic confinement fusion and astrophysics. 

The procedures are vaUd up to high densities. The studies presented in the paper explore 
the density effects in detail within the SCO? picture and show that the density dependencies of 
excited ion populations and of effective rate coefficients cannot be ignored. 

Specific results are presented for light elements up to neon and the computations are 
carried out in an atomic basis of terms (LS-coupled). Such modelling will remain sufficient 
up to about the element argon, beyond which a level basis (intermediate coupling) becomes 
necessary. Heavier elements will be examined in further papers of this series. 

Considerable attention has been given to the generation and assembly of high quality 
fundamental data in support of the SCIR modelhng. Also datasets of fundamental and derived 
data have been specified precisely and codes have been organised following the principles of 
the ADAS Project. The product of the study is the preferred ADAS data for the light element 
ions at this time. 
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